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Abstract 

A wheel or sphere rolling without slipping on the inside of a sphere in a uniform grav- 
itational field can have stable circular orbits that lie wholly above the "equator", while a 
particle sliding freely cannot. 



1 Introduction 

In a recent article |l[ in this Journal, Abramowicz and Szuszkiewicz remarked on an inter- 
esting analogy between orbits above the equator of a "wall of death" and orbits near a black 
hole; namely that the centrifugal force in both cases appears to point towards rather than 
away from the center of an appropriate coordinate system. Here we take a "wall of death" 
to be a hollow sphere on the Earth's surface large enough that a motorcycle can be driven 
on the inside of the sphere. The intriguing question is whether there exist stable orbits for 
the motorcycle that lie entirely above the equator (horizontal great circle) of the sphere. 

In ref. U the authors stated that no such orbits are possible, perhaps recalling the well- 
known result for a particle sliding freely on the inside of a sphere in a uniform gravitational 
field. However, the extra degrees of freedom associated with a rolling wheel (or sphere) 
actually do permit such orbits, in apparent defiance of intuition. In particular, the friction 
associated with the condition of rolling without slipping can in some circumstances have an 
upward component large enough to balance all other downward forces. 

In this paper we examine the character of all circular orbits inside a fixed sphere, for both 
wheels and spheres that roll without slipping. The rolling constraint is velocity dependent 
(non-holonomous), so explicit use of a Lagrangian is not especially effective. Instead we 
follow a vectorial approach as advocated by Milne (Chap. 17) 0. This approach does utilize 
the rolling constraint, a careful choice of coordinates, and the elimination of the constraint 
force from the equations of motion, all of which are implicit in Lagrange's method. The 
vector approach is, of course, a convenient codification of earlier methods in which individual 
components were explicitly written out. Compare with classic works such as those of Lamb 
(Chap. 9) 0, Deimel (Chap. 7) @ and Routh (Chap. 5) §. 

Once the solutions are obtained in sec. 2 for rolling wheels we make a numerical evaluation 
of the magnitude of the acceleration in g's, and of the required coefficient of static friction 
on some representative orbits. The resulting parameters are rather extreme, and the circus 
name "sphere of death" seems apt. 

The stability of steady orbits of wheels is considered in some detail, but completely 
general results are not obtained (because the general motion has four degrees of freedom). 
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All vertical orbits are shown to be stable, as are horizontal orbits around the equator of the 
sphere. We also find that all horizontal orbits away from the poles are stable in the limit of 
small wheels, and conjecture that the a similar condition holds for "death-defying" orbits of 
large wheels above the equator of the sphere. In sec. 4 we lend support to this conjecture 
by comparing to the related case of a sphere rolling within a sphere for which a complete 
stability analysis can be given. 

Discussions of wheels and spheres rolling outside a fixed sphere are given in sees. 3 and 
5, respectively. 

2 Wheel Rolling Inside a Fixed Sphere 
2.1 Generalities 

We consider a wheel of radius a rolling without slipping on a circular orbit on the inner 
surface of a sphere of radius r > a. The analysis is performed in the lab frame, in which the 
sphere is fixed. The z-axis is vertical and upwards with origin at the center of the sphere as 
shown in Fig. |T[ As the wheel rolls on the sphere, the point of contact traces a path that is 
an arc of a circle during any short interval. In steady motion the path forms closed circular 
orbits which are of primary interest here. We therefore introduce a set of axes (x 1 , y', z') that 
are related to the circular motion of the point of contact. If the motion is steady, these axes 
are fixed in the lab frame. 

The normal to the plane of the circular orbit through the center of the sphere (and also 
through the center of the circle) is labeled z' . The angle between axes z and z' is (3 with 
0</3<7r/2. A radius from the center of the sphere to the point of contact of the wheel 
sweeps out a cone of angle 9 about the z' axis, where < < it. The azimuthal angle of the 
point of contact on this cone is called <p, with = defined by the direction of the x'-axis, 
which is along the projection of the z-axis onto the plane of the orbit, as shown in Fig. 0. 
Unit vectors are labeled with a superscript A , so that y' = z' x x' completes the definition 
of the '-coordinate system. 

For a particle sliding freely, the only stationary orbits have (3 = (horizontal circles) or 
(3 = tt/2 (vertical great circles). For wheels and spheres rolling inside a sphere it turns out 
that (3 = or tt/2 also, as we will demonstrate. However, the friction at the point of contact 
in the rolling cases permits orbits with a larger range of 9 than in the sliding case. If (3 = 
or 7r/2 were accepted as an assumption the derivation could be shortened somewhat. 

We also introduce a right-handed coordinate triad of unit vectors (1, 2, 3) related to the 
geometry of the wheel. Axis 1 lies along the symmetry axis of the wheel as shown shown in 
Fig. [I]. Axis 3 is directed from the center of the wheel to the point of contact of the wheel 
with the sphere. The vector from the center of the wheel to the point of contact is then 

a = a3. (1) 

Axis 2 = 3x1 lies in the plane of the wheel, and also in the plane of the orbit (the x'-y' 
plane). The sense of axis 1 is chosen so that the component u\ of the angular velocity vector 
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Figure 1: A wheel of radius a rolls without slipping on a circular orbit inside 
a fixed sphere of radius r. The orbit sweeps out a cone of angle 9 about the 
z'-axis, which axis makes angle (3 to the vertical. The x'-axis is orthogonal 
to the z'-axis in the z-z' plane The angle between the plane of the orbit and 
diameter of the wheel that includes the point of contact with the sphere is 
denoted by a. A right-handed triad of unit vectors, (1, 2, 3), is defined with 
1 along the axis of the wheel and 3 pointing from the center of the wheel to 
the point of contact. 

lj of the wheel about this axis is positive. Consequently, axis 2 points in the direction of the 
velocity of the point of contact, and therefore is parallel to the tangent to the orbit. 

Except for axis 1, these rotating axes are not body axes, but the inertia tensor is diagonal 
with respect to them. We write 

I n = 2kma 2 , J 22 = kma 2 = / 33 , (2) 

which holds for any circularly symmetric disc according to the perpendicular axis theorem; 
k = 1/2 for a wheel of radius a with mass m concentrated at the rim, k = 1/4 for a uniform 
disc, etc. 

The wheel does not necessarily lie in the plane of the orbit. Indeed, it is the freedom 
to "bank" the wheel that makes the "death-defying" orbits possible. The diameter of the 
wheel through the point of contact (i.e., axis 3) makes angle a to the plane of the orbit. 
In general, a wheel can have an arbitrary rotation about the 3-axis, but the wheel will roll 
steadily along a closed circular orbit orbit only if angular velocity component uj% is such that 
the plane of the wheel intersects the plane of the orbit along the tangent to the orbit at the 
point of contact. Hence, for steady motion we will be able to deduce a constraint on u 3 . 
The case of a rolling sphere is distinguished by the absence of this constraint, as considered 
later. 
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Figure 2: The azimuth of the point of contact of the wheel with the sphere 
to the x'-axis is <fi. The unit vector r' is orthogonal to the z'-axis and points 
towards the center of the wheel (or equivalently, towards the point of contact). 
Unit vector 2 = 3 x 1 = z' x r'. 



Since the wheel lies inside the sphere, as shown in Fig. ^, we can readily deduce the 
geometric relation that 

9 - tv + sin _1 (a/r) < a < 9 - sin _1 (a/r). (3) 
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Figure 3: Geometry illustrating the extremes of angle a. 



It is useful to introduce r' = r'r' as the perpendicular vector from the z'-axis to the center 
of the wheel. The magnitude r' is given by 



r' = r sin 6* — a cos a, (4) 
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as shown in Fig. The vector z' x F is in the direction of motion of the point of contact, 
which was defined previously to be direction 2. That is, (f, 2, z') form a right-handed unit 
triad, which is related to the triad (1, 2, 3) by 

z' = — cosal — sin a3, (5) 

and 

r' = 2 x z' = — sin al + cos a3, (6) 

as can be seen from Fig. 1. 

The length r' is negative when the center of the wheel is on the opposite side of the z'-axis 
from the point of contact. This can occur for large enough a/r when the point of contact is 
near the z'-axis, such as when 9^0 and a < or 6 « ix and a > 0. 




Figure 4: Geometry illustrating the vector r cm from the center of the sphere 
to the center of the wheel, and the distance r' — r sin 9 — a cos a from the 
z'-axis to the center of the wheel. 



The force of contact of the sphere on the wheel is labeled F. For the wheel to be in contact 
with the sphere the force F must have a component towards the center of the sphere, which 
will be verified after the motion is obtained. 

The equation of motion of the center of mass of the wheel is 

d 2 r 

where g is the acceleration due to gravity. The equation of motion for the angular momentum 
L cm about the center of mass is 

^ = N cm = ax F. (8) 
at 

We eliminate the unknown force F in eq. (]8|) via eqs. (H) and (0) to find 

l^= 3 3xz + 3x^. (9) 

ma at at z 
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The constraint that the wheel rolls without slipping relates the velocity of the center of 
mass to the angular velocity vector Co of the wheel. In particular, the velocity vanishes for 
that point on the wheel instantaneously in contact with the sphere: 

Vcontact = V cm + UJ X a = 0, (10) 

and hence 



-a3xCo. (IV 



dt 

Multiplying this equation by 3, we find 

Q = _3 x — + to 3 3. (12) 
a 

Equations ©-(0) hold whether the rolling object is a wheel or a sphere. 

The strategy now is to extract as much information as possible about the angular velocity 
Co before confronting the full equation of motion (|9]). The angular velocity can also be written 
in terms of the unit vector 1 along the symmetry axis of the wheel as 

- - dl , , 

u = u l l + lx— . (13) 

This follows on writing Co = uil + Co±, and noting that the rate of change of the body vector 
1 is just dl/dt = Co± x 1, so Coj_ = 1 x dl/dt. Using eq. (fj), the angular momentum can now 
be written as 

L cm — I ■ Co — 2kma 2 u)il + kma 2 l x — . (14) 

dt 

2.2 Steady Motion in a Circle 

To obtain additional relations we restrict our attention to orbits in which the point of contact 
of the wheel with the sphere moves in a closed circle. In such cases the center of mass of the 
wheel (and also the coordinate triad (1, 2, 3)) has angular velocity </> about the z'-axis (and 
no other component), where the dot means differentiation with respect to time. Thus 

v cm = 0z' x r'r = r'02. (15) 

Equation fll~2]) can now be evaluated, yielding 

u=(r'/a)j>i+w 3 3. (16) 

For steady motion there can be no rotation about axis 2; angle a is constant. To find lo 3 we 
now pursue eq. dl3|). 

As argued above, the angular velocity 7 of the triad (1, 2, 3) is 

7 = <pz' = — 0cosal — 0sina3, (17) 

using eq. (|). Then, 

dl - ■ , , 

_ = 7 x 1 = -0sina2 (18) 
dt 
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d2 ■ 

— = 7x2 = 0sinal — 0cosa3 = — 0r', (19) 

(JjV 

and 

d3 

— = 7 x 3 = 0cosct2. (20) 

It immediately follows that 

i x -i = -0sina3. (21) 
dt 

Comparing with eq ( p~3| ) we see that co>3 = — sin a and hence from eq. (|T6| ) we find 

= {r'/a)(j}i - sin a3. (22) 

As anticipated, the rolling constraint specifies how U\ and uj 3 are both related to the angular 
velocity of the wheel about the z'-axis. 

For use in the equation of motion ([|) we can now write 

L = I ■ u = kma 2 [2(r' / a)(j)i - 0sina3], (23) 

and hence, 

— = 2kr'6l — k(fi 2 sin a(2r' + a cos a) 2 — /ca0 sin a3, (24) 

ma at 

using eqs. ( |T8"| - [ZU| ) . Also, by differentiating eq. flip] ) we find 

d 2 r 

Ul A cm 
dt 2 

so that 



r^sinal + r'02 -r'0 2 cosa3, (25) 



d 2 r 

3 x — ^ = -r'01 + r'0 2 sin a2. (26) 
dt 2 

Combining (|9]), fl24"|) and (pBj), the equation of motion reads 

* A d 2 r cm 1 dL 
x 3 — — 3 x 

dt 2 ma dt 

= -(2k + l)r'0i + [(2k + l)r' + fca cosa]0 2 sin a2 + A;a0sina3. (27) 

To evaluate z x 3, we first express z in terms of the triad (f, 2, z'), and then transform to 
triad (1, 2, 3). When the point of contact of the wheel (and hence the r'-axis) has azimuth 
relative to the x' axis, the z axis has azimuth — relative to the f axis. Hence, 

z = sin P cos 0r ' — sin /3 sin 02 ! + cos f3z' (28) 
= —(cos a cos /3 + sin a sin (3 cos 0)1 — sin (3 sin 02 — (sin a cos (3 — cos a sin f3 cos 0)3, 

using eqs. Thus, 

z x 3 = — sin p sin 01 + (cos a cos /3 + sin a sin /3 cos 0)2. (29) 
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The 1, 2 and 3 components of the equation of motion are now 



(2k + l)r r (j) = #sin/3sin0, (30) 
[(2k + l)r' + /cacoso!]0 2 sina = g(cosacos/3 + sin a sin/? cos 0), (31) 

and 

kacf) sin a = 0. (32) 
The cone angle # enters the equations of motion only through r'. 

2.2.1 Vertical Orbits 



From eq. (|32j) we learn that for circular orbits either sin a = or <j) = 0. We first consider 
the simpler case that sin a = 0, which implies that the plane of the wheel lies in the plane 
of the orbit. For a wheel inside the sphere with sin a = 0, we must have a = to satisfy the 
geometric constraint (|3]). Then eq. (|31~D can only be satisfied if cos (3 = 0; i.e., (3 = n/2 and 
the plane of the orbit is vertical. The remaining equation of motion (|30|) now reads 

(2k + l)r'0 = #sm0, (33) 

with r' = r sin 9 — a > 0, which integrates to 

— - — mr (<p — <p Q ) = mgr (1 — cos (p), (34) 

where <f)Q is the angular velocity at the top of the orbit at which = 0. Equation (|34T) 
expresses conservation of energy. The angular velocity u and the angular momentum L cm 
vary in magnitude but are always perpendicular to the plane of the orbit. 

The requirement that the wheel stay in contact with the sphere is that the contact force 
F have component F± that points to the center of the sphere. On combining eqs. (0), (|25|). 
(ID and (H| we find 



2k 

rag sin 02 + m(g cos (b — r'(b 2 )3. (35) 
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The contact force is in the plane of the orbit, so the resulting torque about the center of 
mass of the wheel changes the magnitude but not the direction of the angular momentum. 
On the vertical orbits, axis 2 is tangent to the sphere, and axis 3 makes angle tt/2 — 9 to 
the radius from the center of the sphere to the point of contact. Hence 

F ± = -F 3 sin 9 (36) 

is positive and the orbit is physical so long as the angular velocity 0o & t the peak of the orbit 
obeys 



> 3. (37) 



as readily deduced from elementary considerations as well. 
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The required coefficient fi of static friction is given by \x = F»/F± where 

F,| = ^Ficos6 2 + Fi (38) 

is the component of the contact force parallel to the surface of the sphere. We see that 

H = cot 9s[l + (F 2 /F 3 cos9) 2 , (39) 

which must be greater than cot#, but only much greater if the wheel nearly loses contact 
at the top of the orbit. Hence orbits with tt/A < 9 < tt/2 are consistent with the friction of 
typical rubber wheels, namely fi ^ 1. 

Because a wheel experiences friction at the point of contact, vertical orbits are possible 
with 9 < it/2. This is in contrast to the case of a particle sliding freely on the inside of a 
sphere for which the only vertical orbits are great circles (6 = tt/2). The only restriction 
in the present case is that the wheel fits inside the sphere, i.e., rsin# > a, and that the 
minimum angular velocity satisfy eq. (|37|). 



2.2.2 Horizontal Orbits 

The second class of orbits is defined by <ft — 0, so that the angular velocity is constant, say 
= Q. From eq. ( |30|) we see that sin (3 — and hence (3 = for these orbits, which implies 



that they are horizontal. Then eq. (|3lD gives the relation between the required angular 



velocity f2 and the geometrical parameters of the orbit: 

^9 Q cot a o cot a , 

Q = - = , (40) 

(2k + l)r' + ka cos a (2k + l)r sin 9 — (k + l)a cos a 

recalling eq. (|J). Compare Ex. 3, sec. 244 of Routh |J or sec. 407 of Milne [Q. There are 
no steady horizontal orbits for which a = 0, i.e., for which the wheel lies in the plane of the 
orbit. For such an orbit the angular momentum would be constant, but the torque on the 
wheel would be nonzero in contradiction. 

In the following we will find that horizontal orbits are possible only for < a < tt/2. 

First, the requirement that Q 2 > for real orbits puts various restrictions on the param- 
eters of the problem. We examine these for the four quadrants of angle a. 

1. < a < tt/2. Then cot a > so we must have 

. ka cos a 

r > . (41) 

2k + 1 v ; 

This is satisfied by all r' > and some r' < 0. However, for the wheel to fit inside the 
sphere with < a < tt/2, we can have r' < only for 9 > tt/2 according to eqs. (|f) 
and (|). 

2. tt/2 < a < tt. Then cos a < and cot a < so the numerator of (|40| ) is negative and 
the denominator is positive. Hence Q is imaginary and there are no steady orbits in 
this quadrant. 
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3. —n < a < —tt/2. Then cos a < but cot a > so Q 2 > and r' > and eq. (|40|) 
imposes no to restriction. For the wheel to fit inside the sphere with a in this quadrant 
we must have < 7r/2. 

4. —tt/2 < a < 0. Then cot a < so we must have 

, ka cos a 

r < - a+r < °- (42) 

For the wheel to be inside the sphere with r' < and a in this quadrant we must have 
9 < tt/2. 

To obtain further restrictions on the parameters we examine under what conditions the 
wheel remains in contact with the sphere. The contact force F is deduced from eqs. (0), (f25p 
and ([2D to be 

F/m = (—g cos a + r'VL 2 sin a)l — (g sin a + r'VL 2 cosa)3. (43) 

It is more useful to express F in components along the f and 9 axes where f points away 
from the center of the sphere and 9 points towards increasing 9. The two sets of axes are 
related by a rotation about axis 2: 

1 = -cos(0 - a)? + sin(0 - a)8, 3 = sin(0 - a)r + cos(0 - a)§, (44) 

so that 

F/m = -(r'ti 2 sin0 - gcos0)f - (r'fi 2 cos0 + g sin 0)0 

= -r'Q 2 r + gz. (45) 



The second form of eq. ((45" ) follows directly from elementary considerations. The inward 
component of the contact force, F± = —F r , is positive and the orbits are physical provided 

r'Q 2 > g cot 0. (46) 

There can be no orbits with r' < and < tt/2, which rules out orbits in quadrant 4 of 
a, i.e., for — n/2 > a < 0. 

Using eq. ( ^0|) for Q 2 in eq. ( ^6|) we deduce that contact is maintained for orbits with 
r' > only if 

cot a > [2k + 1 + k(a/r') cos a] cot 0. (47) 

For r' < the sign of the inequality is reversed. 

In the third quadrant of a we have cos a < 0, so inequality ( |4TD can be rewritten with 
the aid of (f|) as 

cota>(l + 2A; r ) cot > cot 0. (48) 

\ 1 + r sin 9/ a |cos a\ J 

However, in this quadrant inequality @ tells us 



cot a < cot 



+ sin" 1 (a/r)l < cot0. (49) 
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Hence there can be no steady orbits with —n < a < —n/2. 

Thus steady horizontal orbits are possible only for < a < n/2. Furthermore, since the 
factor in brackets of inequality (f£?D is roughly 2 for a wheel, this kinematic constraint is 
somewhat stronger than the purely geometric relation (|3|). However, a large class of orbits 
remains with 9 < n/2 as well as 9 > n/2. 

The coefficient of friction \x at the point of contact must be at least F»/F± where F» = \Fq\ 
from eq. (pE5p. (For 9 > n/2 and a near zero the tangential friction Fg can sometimes point 
in the +9 direction.) Hence we need 



> 



\r'Q 2 cos 9 + g sin 6*| 
r'Q 2 sin 9 — g cos 9 



(50) 



The acceleration of the center of mass of the wheel is r'Q 2 , so according to eq. (|l0|) the 
corresponding number of g's is 



cot a 



2k + 1 + k(a/r') cos a 



(51) 



Table [T] lists parameters of several horizontal orbits for a sphere of size as might be found 
in a motorcycle circus. The coefficient of friction of rubber tires is of order one, so orbits 
more than a few degrees above the equator involve very strong accelerations. The head of 
the motorcycle rider is closer to the vertical axis of the sphere than is the center of the wheel, 
so the number of g's experienced by the rider is somewhat less than that given in the Table. 

Figure [| illustrates the allowed values of the tilt angle a as a function of the angle 9 of 
the plane of the orbit, for a/r = 0.1 as in Table |I[ 
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Figure 5: The allowed values of the tilt angle a as a function of the angle 
9 of horizontal orbits for a/r = 0.1. The allowed region is bounded by three 
curves, derived from expressions (|3|), (|40|) and (|16|). 



From eq. ( |^0"D we see that a = tt/2, Q = is a candidate "orbit" in the lower hemisphere. 
On such an "orbit" the wheel is standing vertically at rest, and is not stable against falling 
over. We infer that stability will only occur for Q greater than some minimum value not 
revealed by the analysis thus far. 
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Table 1: Parameters for horizontal circular orbits of a wheel of radius 0.3 
m rolling inside a sphere of radius 3.0 m. The wheel has coefficient k = 1/2 
pertaining to its moment of inertia. The polar angle of the orbit is 9, so orbits 
above the equator of the sphere have 9 < 90°. The plane of the wheel makes 
angle a to the horizontal. The minimum coefficient of friction required to 
support the motion is [i. The magnitude of the horizontal acceleration of the 
center of mass is reported as the No. of g's. 
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2.3 Stability Analysis 

A completely general analysis of the stability of the steady circular orbits found above appears 
to be very difficult. We give a fairly general analysis for vertical orbits, but for horizontal 
orbits we obtain results only for orbits with 9 = tt/2, i.e., orbits about the equator of the 
sphere, and for orbits of "small" wheels. 

We follow the approach of sec. 405 of Milne where it was shown how the steady motion 
of a disk rolling in a straight line on a horizontal plane is stable if the angular velocity is 
great enough. It was also shown that the small oscillatory departures from steady motion 
lead to an oscillatory path of the point of contact of the wheel with the plane. Hence in the 
present case we must consider perturbations that carry the wheel away from the plane of 
the steady orbit. 

The difficulty is that there are in general four degrees of freedom for departures from 
steady motion: the axis of the wheel can be perturbed in two directions and the angular 
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velocity can be perturbed as well as the angle 9 to the point of contact. However, the 
procedure to eliminate the unknown force of contact from the six equations of motion of a 
rigid body leaves only three equations of motion. We will obtain solutions to the perturbed 
equations of motions only in special cases where there are in effect just two or three degrees 
of freedom. A more general analysis might be possible using the contact force found in 
steady motion as a first approximation to the contact force in perturbed motion, but we do 
not pursue this here. 

A wheel rolling with a steady circular orbit on a plane can suffer only three types of 
perturbations and the results of an analysis are reported in Ex. 3, sec. 244 of Routh ||. For 
a sphere rolling within a fixed sphere the direction of what we call axis 3 always points to the 
center of the fixed sphere so there are only two perturbations to consider and the solution is 
relatively straightforward, as reviewed in sec. 4 below. The stability of horizontal orbits of 
rolling spheres lends confidence that stable orbits also exist for wheels. 

2.3.1 Vertical Orbits 

We define the (F, y', z') coordinate system to have the F-axis vertical: x' = z. In steady 
motion we have 

a = 0, 1 = — z', and 3 = F = x'cos0 + y'sin0, (52) 

where is the azimuth of the center of the wheel from the x'-axis. Thus = at the top of 
the orbit. To discuss departures from steady motion in which the 1-axis is no longer parallel 
to the z'-axis, it is useful to have a unit triad (F, s', z') defined by eq. ( 152] ) and 

s' = z' x f ' = — x' sin + y' cos 0, (53) 

with defined as before. See Fig. |6]. The surface of the sphere at the point of contact is 
parallel to the s'-z' plane. Axes F and z' rotate about the z'-axis with angular velocity 0, 
so that 

- = 0s, and ^ = -0r. (54) 

The perturbed 1-axis can then be written 

1 = e r F + e s s' — z', with \e r \ , \e s \ <C 1, (55) 

where throughout the stability analysis we ignore second-order terms. Writing 

3 = F + 5 8 s' + 5 z z', with \S S \ , \S Z \ < 1, (56) 

the condition 1-3 = requires that S z = e r . Then to first order, 

2 = 3x1 = -S s r' + s' + e s z. (57) 

We expect that vector 2 will remain parallel to the surface of the sphere even for large 
departure from steady motion, so 2 must remain in the s'-z' plane. Hence, 5 S = 0, and 

3 = F + e r z'. (58) 
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Also, we can identify a as the tilt angle of the 3-axis to the r'-s' plane, so that 

a = e r . (59) 

The analysis proceeds along the lines of sec. 2.1 except that now we express all vectors 
in terms of the triad (f, s', z'). To the first approximation the angular velocity of the wheel 
about the 1-axis is still given by U\ = (r'/a)<p. From eqs. ( |5*^ ) and flSUP we find 

^ = (e r - e s 4>)v' + (e r + e s )s', (60) 
1 x ^ = (e r (p + e s )i'- (e r - e s 0)s', (61) 

so that eq. (O) yields 

£ = cud + i x ^ = [(1 + r'/a)e r + e s ]r' - [e r - (1 + r'/a)e s 0]s' - (r'/a)<j>z'. (62) 
Then eq. ([lj|) tells us 

— = 2kaoo 1 l + kal x ^ = fc[(2r' + a)e r + ae s ]f - k[ae r - (2r + a)e s 0]s' - 2kr'<j)z', (63) 
so that to first order of smallness 

— -^ = k[2(r' + a)e r <p + (2r + a)(e r (f) - e s <p 2 ) + ae s ]r 

TYIjCL Cub 

-k[ae r - (2r' + a)(e r 2 - e s '(p) - 2{r' + a)e s 0]s' 

-2k(r'4> + r^)z'. (64) 



11 



In this we have noted from eq. (f4|) that f ' = r6 sin 9 to first order, and that 6 is small. Next, 
dr rr 



dt 

Then to first order, 
d 2 r 



a3 x Co ~ av' x Co — r'0s' — [ae r — (r + a)e s (p}z. (65) 



r'0 2 r' + ( r '0 + +'(/))&' - [ae r - (r' + a)(e s + e s 0)]z', (66) 



dt 2 
so that 



3 x — ^ = f ' + e r z' x 



cm 



dt 2 v T ' dt 2 

-r'e r 0r' + [ae r - r'e r 2 - (r' + a)(e s <j) + e a <j>)]s' + (r'0 + f'0)z' (67) 

Also, 

3 x z = (f' + e r z') x (cos <fii' — sin <ps') = e r sin 0r' + e s cos 0s' — sin <pz' . (68) 
The r', s' and z' components of the equation of motion @ are then 

= [(2k + l)r' + ka\e r (f) + 2k(r' + a)e r (p — ge r sin <p — k(2r' + a)e s <fi 2 + kae s , (69) 

= [{2k + l)r + ka}e r (p 2 — {k + l)ae r — ge r cos^> 

+ (2/e + l)(r' + a)e s + [(2k + l)r' + (/c + l)a]e s 0, (70) 

= (2k + l)(r'4> + f'<f)) -gsmcf). (71) 

If the perturbations e r , e s and r' are set to zero eqs. (|69|) and (|70|) become trivial while 
eq. (|7TD becomes the steady equation of motion (p3|) . 

The general difficulty with this analysis is that there are only three equations, ( |69| - |7T| ), 
while there are four perturbations, e r , e s , <p and 9. The perturbation 9 appears only in 
eq. (|7l|) via r'; its effect on r' leads only to second-order terms in eqs. (|69|-[f0|). If we could 



neglect the terms in in eqs. (p9| - [70[) then these two equations would describe only the 
perturbations e r and e s to first order and a solution could be completed. 



Therefore we restrict our attention to the top of the orbit, = 0, where eq. (|7lD tells 
us that = to leading order. The angular velocity 0o at this point is a minimum so the 
gyroscopic stability of the wheel is the least here. Hence if the orbit is stable at = it will 
be stable at all 6. 



The forms of eqs. ( |69"D and (|70| ) for = indicate that if e r and e s are oscillatory then 
they are 90° out of phase. Therefore we seek solutions 

e r = e r cos tut, e s = e s sin cut, (72) 

where lo now represents the oscillation frequency. The coupled equations of motion then 
yield the simultaneous linear equations 

2k(r' + a)(j)QUJ€ r + [kau 2 + k(2r' + a)(fil}e s = 
{(k + l)au 2 + [(2k + l)r' + ka]<p 2 - g}e r + (2k + l)(r' + a)j) ue a = 0. (73) 
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These equations are consistent only if the determinant of the coefficient matrix vanishes, 
which leads to the quadratic equation 

Au 4 - Bu 2 - C = 0, (74) 

with solutions 



B ± VB 2 + 4AC 



where 



and 



" = —^A • < 75 > 

A = k(k + l)a 2 , (76) 
B = kag + k[(2k + l)(2r /2 + a 2 ) + (4k + l)ar> 2 , (77) 



C = k(2r' + a) {[{2k + l)r' + fca]$ - g) 2 . (78) 

Since A and B are positive there are real, positive roots whenever B 2 + 4AC is positive, 
i.e., for C > — B 2 /4A. In particular, this is satisfied for positive C, or equivalently for 

* > (2t + iy + fa - (79) 

However, this is less restrictive than the elementary result ( PTp that the wheel stay in contact 
with the sphere! All vertical orbits for which the wheel remains in contact with the sphere 
are stable against small perturbations. 



The stability analysis yields the formal result that if (0, 0) = (0, 0) then u> — J g/(k + l)a. 
We recognize this as the frequency of oscillation of a simple pendulum formed by suspending 
the wheel from a point on its rim, the motion being perpendicular to the plane of the wheel. 



2.3.2 Horizontal Orbits 

We expect the stability analysis of horizontal orbits to be nontrivial since we have identified 
steady orbits that are "obviously" unstable. 

The spirit of the analysis has been set forth in the preceding sections. For horizontal 
orbits the (x 1 , y', z') coordinate system can be taken as identical with the (x, y, z) system, so 
we drop symbol ' in this section. We introduce a triad (r, s, z) with f being the perpendicular 
unit vector from the z-axis toward the center of the wheel. Then s points in the direction of 
the motion of the center of the wheel in case of steady motion. 

It is also useful to introduce a unit triad that points along the (1, 2, 3) axes for steady 
motion. The s axis already points along the 2 axis for steady motion, so we only need define 
t as being along the direction of 3, and u as being along the direction of 1 for steady motion, 
as shown in Fig. [7]. Then, (s, t, u) form a right-handed unit triad. The vertical, z, is then 
related by 

z = — sin aot — cos ao u , (80) 

where ao is the angle of inclination of the wheel to the horizontal in steady motion. The 
triad (s, t, u) rotates about the z-axis with angular velocity <p, so that 

ds ■ 

— = <pz x s = —(f) cos a t + 4> sin a u ; (81) 
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and 



dt ■ 

— = <pcosaos, 
dt 



du ■ 

— = -ipsinaos. 
dt 



(82) 
(83) 




Figure 7: For horizontal orbits of a wheel rolling inside a sphere the (x, y, z) 
axes are identical with the (x', y', z') axes. The r-s plane is horizontal. The 
axes u, s and t are along the unperturbed directions of the 1, 2 and 3 axes, 
respectively. Axes t and u lie in the vertical plane r-z. 



We now consider small departures from steady motion. The 1-axis deviates slightly from 
the u-axis according to 

i = e s s + e t t + u, |e a |,|e t |<l. (84) 

The 3-axis departs slightly from the t-axis, but to the first approximation it remains in a 
vertical plane, i.e., the t-u plane. Then we have 

2 = s — e s u, and 3 = t — e t u. (85) 

With the above definitions the signs of angles a and e t are opposite: 

Act = —e t , a = —e t . (86) 

To first approximation the component u\ of the angular velocity of the wheel about its 
axis remains uj\ = (r'/a)<fi. Then 

— = (— <j) sin a + e s — e t 4> cos a )s — (e s <j) cos a — e t )i + e s 4> sin a u, (87) 

(JjL 

1 x — = (e s <j) cos a — e t )s — (0 sin cto + e s — e t 4> cos ao)i + e t <j) sin aou, (88) 
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so that 



UJ = CJxl + lX — 

dt 



[(r'/a + cos a )e s <fi - e f ]s 

— [0sin«o — e s — (r'/a + cos«o)e t 0]t 

+ (r'/a + e t sin « o )0u, (89) 



and 



L , s , s dl 

= 2ka<jj\ \ + real x — 

ma at 

= fc[(2r' + a cosao)e s (j) — ae t ]s 

—k[a<p sin «o — ae s — (2r' + a cos «o)e t 0]t 

+A;(2r' + ae t sina o )0u. (90) 

Then to the first approximation 

— = — k[(2r' + a cos «o)0 2 sin a — 2(r' + a cos «o)e s 0) 

ma at 

— (2r' cos a + a cos 2a o )e t 2 + ae t ]s 
—k[a(j) sin a + (2r' + a cos a o )e s 2 cos a — °£ s — 2(r' + a cos a )e t (f>)]t 
+k[2r r (j) + 2f'0 + (2r' + a cos a )e s (j) 2 sin a ]u. (91) 

Unlike the case of vertical orbits, for horizontal orbits the factor <p has no zeroeth-order 
component and we neglect terms like e<p. 
Similarly 

— — = a3 x a; = a(t — e t u) x to 

(JjL 

= r'(f>s — [(r' + a cos a )e s (f) — ae t ]u, (92) 
d 2 r cm ... 

— = [r'0 + r'0 + (r' + a cosa o )e s 2 sina — ae t 4>sma }s 
at 2 

—r'4> 2 cos otot + [r'0 2 sin «o — (r' + a cos ao)e s + ae't]u, (93) 



and 



a d r cm /? /. \ d r cm 

3x ^ = (t 7 etu)x ^ 

= [r'</> 2 sin a — (r' + a cos ao)e s _ r'e t <fi 2 cos «o + a Q]s 

— [r'0 + r'0 + (r' + acosa)e s 2 sina o — aet0 sin ao]u. (94) 

We also need 

3 x z = (t — e 4 u) x (— sin aot — cos aou) = — (cosao + Q sin ao)s. (95) 
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The s, and t and u components of the equation of motion (^|) are 

= [(2k + l)r' + /ca cos ao]0 2 sin ao — gcosa — (2k + l)(r' + a cos ao]e s (j) 

— [(2k + l)r'cosa + ka cos 2ao\e t cf) 2 — ge t suiaQ + (k + l)ae't, (96) 

= kacf) sin cio + k(2r' + a cos ao)e s (j) 2 sin a — kae s — 2k(r' + a cos ao)et<f>, (97) 

and 

= (2k + l)(r> + r>) + [(2k + l)r' + (fc + l)acosa o ]es0 2 sina o - ae t <psma . (98) 

The leading terms of these three equations are just eqs. (|30|)-(|32|) for j3 = 0. Therefore we 
can write <fi = Q + 5 where Q is the angular velocity of the steady horizontal orbit and 5 is 
a small correction. 

Although the derivative of r', 

r' = r6 cos 6 + aa sin a — r ^ cos — a ^t sin ao, (99) 

appears only in eq. (^1), in general the perturbation 9 is not decoupled from e s and et as 
was the case for vertical orbits. Thus far, we have found a way to proceed only in somewhat 
special cases in which the 9 perturbation can be ignored, as described in sees. 2.3.3 and 2.3.4. 



2.3.3 Orbits Near the Equator 

It appears possible to carry the analysis forward for the special case 9 = ir/2, the orbit on 
the equator of the sphere. This case is, however, of interest. 

Assuming 9q = n/2 the equations of motion (|96|-|98|) then provide three relations among 
the three perturbations 5, e s and e t . For this we consider only the first-order terms, noting 
that <p 2 w tt 2 + 2Q5 and 



r 



' = r sin 9 — a cos a ~ r' + r A9 cos 9q + aAa sin «o = r' — ae t sin «o, (100) 



where r = r — acosao for 9 = n/2, recalling eq. (p^). Also, from the form of eqs. (|96|-|9"8|) 
we infer that if the perturbations are oscillatory then 5 and e s have the same phase which is 
90° from that of e t . Therefore we seek solutions of the form 

5 = 5 sin ut, e s = e s smujt, and e t = e t cos out, (101) 

where io is the frequency of oscillation. The first-order terms of the differential equations 
(|96| - |98D then yield the algebraic relations 

= 2f2sinao[(2/c + l)r' Q + kacosao\u5 — Q(2k + l)(r' + a cos ao)u)e s 

— {VL 2 [(2k + l)r' cos«o + (k + sin 2 ao)a] — g sina + (k + l)au; 2 }e t , 
= —ka sin a u 2 5 + [kVL 2 cos a (2r' + a cos a ) + kau 2 ]e s + 2kQ(r' + a cos a )u)e t , (102) 
= -(2k + l)r' LU 2 5 + n 2 sin a [(2k + l)r' + (k + l)a cos a ]e s + 2(k + sin a aue t . 
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These equations have the form 



A n w5 + A 12 ue s + (A 13 + B 13 u 2 )e t = 

A 21 u 2 5 + {A 22 + B 22 co 2 )e s + A 23 toe t = (103) 

A 31 u 2 S + A 32 e s + A 33 ue t = 

To have consistency the determinant of the coefficient matrix must vanish, which leads 
quickly to the quadratic equation 

Aoo 4 - Boo 2 - C = 0, (104) 

where 

A = B 13 B 22 A 31l (105) 
B = A n B 22 A 33 + A l2 A 23 A 31 + B 13 A 21 A 32 - A 13 B 22 A 31 - B 13 A 22 A 3l - A 12 A 2l A 33l (106) 
and 

C = A n A 22 A 33 + A l3 A 21 A 32 - A 13 A 22 A 3l - A U A 23 A 32 . (107) 

From numerical evaluation it appears that A, B and C are all positive for angular ve- 
locities Q that obey eq. (0). That is, all steady orbits at the equator of the sphere are 
stable. There is both a fast and slow oscillation about steady motion for these orbits, an 
effect familiar from nutations of a symmetric top. 



2.3.4 Small Wheel Inside a Large Sphere 

The analysis can also be carried further in the approximation that the radius a of the wheel 
is much less than the radius r of the fixed sphere. In this case the perturbation in angle 9 
of the orbit is of higher order than the perturbations in azimuth and in the angles e s and 
e t related to the axes of the wheel. A solution describing the three first-order perturbations 
can then be obtained. 

For the greatest simplification we also require that 

a r' w r sin 8 Q . (108) 

Thus we restrict our attention to orbits significantly different from the special cases of motion 
near the poles of the fixed sphere. 

In the present approximation the first-order terms of the perturbed equations of motion 
(H^D are 

2VL5 sin a = fie, + (n 2 cos a + ^k + l) r ' ) Q ' (1 ° 9) 

(HO) 



Q sin a 
and 

5 = — e s fl 2 sin a . (Ill 
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Inserting ( FTP ) into ( 111 ) we can integrate the latter to find 

5 = -ne t . 

Using this and the derivative of ( |110[ ) in ( |109| ) we find that e t obeys 



Q 2 sin a (cos a + 2 sin a ) 



gsm ao 
(2fc + l)r 



e t = 0. 



(112) 



;ii3) 



The the frequency u of the perturbations is given by 



VL 2 sin ao(cos ao + 2 sin ao) + 



flfsin a o2 



(2k + iy 



fi 2 tana (l + sin2a ), (114) 



using eqs. (|40|) and (|10|). 

Thus all orbits for small wheels are stable if condition ( 108|) holds. We conjecture that 
orbits for large wheels are also stable if (|108|) is satisfied. 

For steady orbits that lie very near the poles, i.e., those that have r' Q < a, we conjecture 
that the motion is stable only for Q greater than some minimum value. For a wheel spinning 
about its axis on a horizontal plane the stability condition is 



Q 2 > 



(2fc + l)a 



(115) 



See, for example, sec. 55 of Deimel H. However, we have been unable to deduce the gener- 
alization of this constraint to include the dependence on r and 9 for small r sin 6*0. 



3 Wheel Rolling Outside a Fixed Sphere 

Equations (fl|)- (|32]) hold for a wheel rolling outside a sphere as well as inside when the 
geometric relation @ is rewritten as 

6<a<n + 6. (116) 

We expect no vertical orbits as the wheel will lose contact with the sphere at some point. 
To verify this, note that the condition sin a = (from eq. fl32|)) implies that a = tt when 
the wheel is outside the sphere. Then eqs. ( |34| - |36| ) indicate, for example, that if the wheel 
starts from rest at the top of the sphere it loses contact with the sphere when 

cos0 = 2 - . (117) 

The result for a particle sliding on a sphere (k = 0) is well known. 

For horizontal orbits, eqs. (^0| - ^5D are still valid, but the condition that friction have an 



outward component is now 

r'ft 2 <£cotfl, (118) 
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and hence 

cot a < (2k + 1 + k(a/r') cos a) cot 9. (119) 

Equation fl40|) can be satisfied for a < n/2 so long at the radius of the wheel is small enough 
that (2k + l)r' + ka cos a is positive. We must have 9 < w/2 to have a < it/2 since a > 9, so 
horizontal orbits exist on the upper hemisphere. A particular solution is a = ir/2 for which 
Q = 0; this is clearly unstable. 

There is a class of orbits with 9 < n/2 and a very near it + 9 that satisfy both eqs. ([EJ) 
and ( |119| ). These also appear to be unstable. 

The stability analysis of the preceding section holds formally for wheels outside spheres, 
but the restriction there to the case of 9 = 90° provides no insight into the present case. 



4 Sphere Rolling Inside a Fixed Sphere 

The case of a sphere rolling on horizontal orbits inside a fixed sphere has been treated by 
Milne @. For completeness, we give an analysis for orbits of arbitrary inclination to compare 
and contrast with the case of a wheel. 

Again the axis normal to the orbit is called z', which makes angle (3 to the vertical z. The 
polar angle of the orbit about z' is 9 and is the azimuth of the point of contact between 
the two spheres. The radius of the fixed sphere is r. 

The diameter of the rolling sphere that passes through the point of contact must always 
be normal to the fixed sphere. That is, the "bank" angle of the rolling sphere is always 
9 — 7i / 2 with respect to the plane of the orbit. 

The rolling sphere has radius a, mass m and moment of inertia kmc? about any diameter. 
The angular momentum is, of course, 

L cm = kma 2 uj, (120) 

where u is the angular velocity of the rolling sphere. 

We again introduce a right-handed triad of unit vectors (1, 2, 3) centered on the rolling 
sphere. For consistency with the notation used for the wheel, axis 3 is directed towards the 
point of contact, axis 2 is parallel to the plane of the orbit, and axis 1 is in the 3-z' plane, 
as shown in Fig. [|[ In general, none of these vectors are body vectors for the rolling sphere. 
The center of mass of the rolling sphere lies on the line joining the center of the fixed sphere 
to the point of contact, and so 

r cm = (r - a)3 = r'3, (121) 

Equations (^-pD that govern the motion and describe the rolling constraint hold for the 
sphere as well as the wheel. Using eqs. ( |120| ) and (|121|) we can write eq. ([|) as 



bf = 9 3xz + r'3x5. (122) 
dt dt 2 K ' 

We seek an additional expression for the angular velocity uj of the rolling sphere, but we 
cannot use eq. (|T3|) since we have not identified a body axis in the sphere. However, with 
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Figure 8: Geometry illustrating the case of a sphere rolling without slipping 
on a circular orbit perpendicular to the z'-axis inside a fixed sphere. The 3- 
axis is along the line of centers of the two spheres, and passes through the 
point of contact. The 2-axis lies in the plane of the orbit along the direction 
of motion of the center of the rolling sphere, and axis 1 = 2 x 3 is in the 3-z' 
plane. 



eq. ( |121[ ) the rolling constraint ( |T2"D can be written 



r' d3 

u = 3x— + u 3 3. (123) 

a at 

We can now see that 003 = uj ■ 3 is a constant by noting that 3 • duj/dt = from eq. (|122|) . 
and also uj ■ d3/dt = from eq. (|123|) . The freedom to chose the constant angular velocity 



UJ3 for a rolling sphere permits stable orbits above the equator of the fixed sphere, just as 
the freedom to adjust the bank angle a allows such orbits for a wheel. 
Taking the derivative of eq. ( |123j ) we find 

du r'~ d 2 3 d3 

A=^ 3X dF + Ws A- (124) 
so the equation of motion ( |122j ) can be written 

(k + l)r'3 x 44 - kawj^- =jzx3. (125) 
at 2 dt 

Milne notes that this equation is identical to that for a symmetric top with one point fixed 
M, and so the usual extensive analysis of nutations about the stable orbits follows if desired. 
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We again restrict ourselves to circular orbits, for which the angular velocity of the center 
of mass, and of 1, 2 and 3 is 0z' where the z'-axis is fixed. Then with 



we have 



z' = -sin 01 + cos 03, (126) 



^ = 0z' x 3 = sin 02, (127) 



dt 2 
and hence, 



d 2 3 ..... ..... 

2 sin 0z' x 2 + sin 02 = -0 2 sin cos 01 + sin 92 + 2 sin 2 03, (128) 



d 2 3 

3 x — = -0 sin 01 -0 2 sin0cos02. (129) 



With these the equation of motion ( |125| ) reads 

(k + l)r'0sin0i + [(A; + l)r'0 2 cos0 + A;ac;30]sin02 = -gz x 3. (130) 
We can use eq. fl29|) for z x 3 if we substitute a = 9 — it /2 for the rolling sphere: 

z x 3 = — sin (3 sin 01 + (sin cos /3 — cos sin /? cos 0)2. (131) 
The components of the equation of motion are then 

(k + l)r'0 sin = sin /3 sin 0, (132) 

[(k + l)r'0 2 cos + kauj 5 (p] sin = g cos sin /3 cos — g sin cos /3. (133) 

The two equations of motion are not consistent in general. To see this, take the derivative 
of eq. (|133|) and substitute from eq. (|132|) : 

kauj^ sin (3 sin = —3(k + l)r'0 cos sin f3 sin 0. (134) 

While this is certainly true for /3 — (horizontal orbits), for nonzero /3 we must have 0cos0 
constant since 0J3 is constant. Equation (|134|) is satisfied for 9 = tt/2 (great circles), but for 



arbitrary we would need constant which is inconsistent with eq. (|132|). Further, on a 



great circle eq. ( |133| ) becomes kau^cp = —g cos (3. This is inconsistent with eq. ( |132| ) unless 
(3 = tt/2 (vertical great circles) and ui% = 0. 

In summary, the only possible closed orbits for a sphere rolling within a fixed sphere are 
horizontal circles and vertical great circles. 

We remark further only on the horizontal orbits. For these = Q is constant according 
to eq. ( |132|) . Equation ( |133|) then yields a quadratic equation for Q: 

(k + l)r'Q 2 cos + kau 3 Q + g = 0, (135) 

so that there are orbits with real values of Q provided 

(kau 3 ) 2 >4(k + l)gr' cos 9. (136) 



21 



This is satisfied for orbits below the equator (9 > tt/2) for any value of the "spin" 003 of the 
sphere (including zero), but places a lower limit on \us\ for orbits above the equator. For the 
orbit on the equator we must have Q = —g/{kaui 3 ) so a nonzero o> 3 is required here as well. 
The contact force F is given by 

F/m = (g + r'Q 2 cos 9) sin 91 - (r'Q 2 sin 2 9 - g cos#)3, (137) 



using eqs. ([7p and ( |134 ). For the rolling sphere to remain in contact with the fixed sphere 



there must be a positive component of F pointing toward the center of the fixed sphere. 
Since axis 3 is radial outward from the fixed sphere, we require that F 3 be negative, and 
hence 

r'Q 2 sm 2 9 > gcos9. (138) 

This is always satisfied for orbits below the equator. For orbits well above the equator this 
requires a larger value of |cc^3 1 than does eq. ( |136| ). To see this, suppose u 3 is exactly at the 
minimum value allowed by eq. ( |136| ), which implies that Q = — kau 3 / {2{k + l)r'cos0). Then 
eq. ( |138[ ) requires that tan 2 9 > k + 1. So for k = 2/5 and at angles 9 < 50° larger values of 
|co"3 j are needed to satisfy eq. (|136|) than to satisfy eq. ( |136|) . However, there are horizontal 



orbits at any 9 > for \lu 3 \ large enough. 

5 Sphere Rolling Outside a Fixed Sphere 

This case has also been treated by Milne ||. A popular example is spinning a basketball on 
one's fingertip. 

Equations eq. ( [L35|) and ( |136| ) hold with the substitution that r' = r + a. The condition 
on the contact force becomes 

r'Q 2 sm 2 9 < gcos9, (139) 

which can only be satisfied for 9 < ir/2. While eq. (|136f) requires a large spin 1 0^3 1 , if it is 
too large eq. ( 139 ) can no longer be satisfied in view of the relation ( |135| ). For any case in 



which the orbit exists a perturbation analysis shows that the motion is stable against small 
nutations 0. 
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